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= I. INTRODUCTION 


A. BACKCROUND 


Despite technological advances in ship silencing, the 
acoustic Signature produced by a submarine's main propulsion 
system and auxiliary machinery remains the single nost 
easily exploited tool for detecting, tracking and localizing 
the submarine. The transmission of sound through the ocean 
is predictable yiven a sufficiently well-defined knowledge 
of the environment, the source and the receiver. The trans- 
misSion of sound in the sea depends on the sound velocity 
profile (SVP). Sound velocity in turn depends on the 
temperature, salinity and depth of the water. Transmission 
phenomena can become extremely complex, even for relatively 
Simple sound velocity profiles. 

Acoustic transmission loss models have been developed 
for implementation on digital computers. The models gener- 
ally take one of two approaches: either solving the 
acoustic wave egquaticn or employing ray-tracing technigues. 


Both approaches require assumptions, sSimplificaticns and 


approximations. For real situations some models will 
predict transmission loss. more accurately than others 
depending on the ccmputational technigue involved, the 


method of solution and the extent to which the model accepts 
real conditions for Loundaries at the sea surface, the sea 
floor and between different water masses. 

The situation of interest here is one of the nore 
complicated cases: the shallow sound channel. Shas. s2tu— 
ation develops in higher latitudes of the world's oceans, 
where precipitation is heavy and evaporation is small. A 


sound velocity profile reveals the presence of a mixed layer 


resulting from winds and turbulence which may extend from 
the surface to depths of several tens of meters. Below the 
mixed layer a shallcw sound channel (SSC) is formed bya 
combination of temperature and salinity gradients, resulting 
in a sound speed minimum (channel axis) at roughly 80 to 200 
meters. The lower edge of the SSC is bounded by an area of 
rapidly increasing scund speed. Below the SSC the combined 
effects of decreasing temperature and increasing fressure 
can also result in a second or deep sound channel (DSC) with 
an axis at about 800 to 1000 meters in mid-latitudes. Sound 
channels act aS acoustic waveguides, trapping certain 
frequencies and allowing them to propagate especially well, 
while other frequencies "leak" out of the channel and propa- 
gate poorly. A more complete discussion of the theoretical 
background for this paper is contained in Chapter II. 

The various layers and channels discussed above occupy 
the tcp several hundred to one thousand meters of the 
oceans. Submarines operate within these depths. These 
features of the sound velocity profile strongly affect the 
propagation of the sukmarine's acoustic Signature, and the 
active and passive operation of the submarine's sonar equip- 
ment and that of its hunters. Antisubmarine warfare (ASW) 
operators must understand these effects and have some means 
of predicting them. With reliable predictions oferators 
will be able to select frequencies and modes of operaticn to 
maximize the effectiveness of their sonar equipment in 


detecting and tracking the submarine. 


Bek OBJECTIVE 


This study seeks to use two acoustic transmission loss 
models to predict the behavior of sound in typical shallow 
sound channels. Both models are available at the Naval 
Postgraduate School (NPS), accessed through the IBM 3033 


computer. 
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The first model, the split-step Parabolic Equation (PE) 
model, uses a parakolic approximation to the basic wave 
equation and solves for pressure levelS asa function of 
range. It is considered to be a low-frequency (below 800 
Hz) mcdel, since higher frequencies reyuire large amounts of 
computational time. The PE model, a full-wave solution, is 
considered a benchmark against which other transmission loss 
models are tested. It is not commonly used in suppcrt of 
fleet operations due to the large computational times 
required. 

The second model, the Fast Asymptotic Coherent 
Transmission loss model version 9H (FACT 9H), employs ray- 
tracing technigues to predict the amount of acoustic energy 
acriving at a receiver as a function of range. This nodel 
is less rigorous in predicting the effects of different 
types of boundaries and it does not allow for variations of 
the sound velocity profile with range. It is simpler and 
reguires less computational time. FACT is generally consid- 
ered more appropriate than the PE model for freguencies 
above 800 Hz. The FACT model is widely used in acoustic 
transmission loss predictions provided to fleet users in the 
ea. NAVY. Both the PE model and FACT 9H are available at 
the Naval Postgraduate School, accessed through the IBM 3033 
computer. 

In 1984 Dosso and Chapman [Ref. 1] reported an experi- 
ment on propagation loss that they had conducted in the 
northeast Pacific Ocean. Sound velocity profiles indicated 
the presence of a shallow scund channel overlying a deep 
sound channel. They placed sources of known magnitude in 
the shallow sound channel, with receivers in the shallow 
sound channel and the deep sound channel. They measured 
propagation losses over ranges extending from zero to 35 km 
for propagation both within the SSC and across the Loundary 


with the DSC. Using a split-step parabolic equation nedel, 
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Dosso and Chapman were able to obtain good agreement between 
their experimental results and the predictions of the model. 
A full discussion of Dosso and Chapman's experiment is 
contained in Chapters trer 

This study examines the Dosso and Chapman experiment, 
then attempts to reproduce their results using the Naval 
Postgraduate School's PE and FACT 9H models. Chapter IV 
compares these predictions with the Dosso and Chapman 
results, supporting the use of the two models to predict the 
optimum freyuency for sound propagation within a_ shallow 
sound channel. The models are then applied to two sound 
velocity profiles of interest to the Naval Postgraduate 
School's Environmental Acoustic Research Group. These 
profiles were obtained under summer conditions ina deep 
ocean basin (4000 meters) at high latitude. Chapter V 
describes the results of the NPS PE model as applied to 
these prefiles. Chapter VI discusses the FACT 9H model and 
its predictions for the Dosso and Chapman profiles and the 
high-latitude summer frofiles. Chapter VII then summarizes 
this study's results and conclusions and presents recommen- 


dations for further investigation. 
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A. SCUND IN THE SEA 


The manner in which sound is transmitted through a 
medium depends on the velocity with which it propagates. 
Because of the greater density of seawater, sound is trans~ 
mitted much faster inthe ocean than in the air. The 
velocity of sound at any given point in the ocean is deter- 
mined by the temperature, Salinity and pressure. Since 
pressure varies linearly with water depth, these two quanti- 
ties can be treated as interchangeable. Temperature and 
depth are generally the dominant factors affecting the speed 
of sound. Thus the velocity of sound in the ocean varies in 
ali dimensions: herizontally, vertically and temporally 
[Ref. 2]. 

The sound velocity can be measured directly using a 
velocimeter or indirectly using a bathythermograph (8T) to 
measure temperature as a function of depth. Numerous 
researchers have provided empirical formulas for calculating 
velocity from temperature, salinity and depth [Ref. 3]. In 
either case the result 1s a record of the sound velocity as 
it changes with depth, or sound velocity profile (SVP). 

Snell's Law of acoustics is exactly analogous to Snell's 
Law of optics. When a ray of light or sound passes from one 
medium intc another, the ray is bent or refracted toward the 
medium in which it travels at the slower speed. The ocean 
may be thought of as composed of many layers of water, in 
each of which sound has a different velocity. Where the SVP 
Shows a velocity minimum, sound is refracted toward that 
depth, resulting in a focusing of energy at that velocity 
minimum. This depth is often referred to as a sound channel 


axis. 
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A typical sound velocity profile like that of Figure 2.1 
[Ref. 3] shows a surface layer of Seawater which is subject 
to the actions of sun, wind and precipitation. This layer 
is often well-mixed, but its temperature, salinity and 


thickness may vary greatly from place to place and from day 


to day. Typically the surface layer shows a constant or 
Slightly increasing sound velocity profile. The surface 
layer may extend to a depth of several tens of meters. At 


some foint below tne surface layer a strong thermal gradient 


may show up. Temperature and therefore sound velccity 
decrease Sharply with depth, creating a seasonal thermo- 
efine:. Deeper still lies the main thermocline, formed by 


another strong negative thermal gradient. Below this is the 
deep isothermal layer, which extends to the sea floor. rn 
this layer salinity is essentially constant, as is tempera- 
ture. Pressure and velocity increase very nearly linearly 
with depth. The layers bounded above by the main thermo- 
cline and below by the increasing velocity gradient are 
Called the deep sound channel (DSC). Sound can be trapped 
in this channel and refracted repeatedly toward the velocity 
Hinimum, which is called the DSC axis. This axis may occur 
at depths ranging frcem several hundred to over 1000 meters. 

Certain environmental conditions can cause a sound 
velocity profile like that shown in Figure 2.2. A combina- 
tion of winds, high precipitation and low evaporation such 
aS 1s found in the northeastern Pacific Ocean can lead to 
the formation of a secondary or shallow sound channel (SSC) 
just kelow the seasonal thermocline. Just as sound can be 
trapped in the DSC, sound can also be trapped in the SSC 
fRef. 1]. 
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Be. ACOUSTIC WAVEGUIDES 


Sound channels act as acoustic waveguides. The ray 
theory cf sound propagation predicts that sound can be 
refracted and trapped in ducts. Trapping depends’ on the 
thickness of the duct and the wavelength of the sound; if 
the wavelength is "too large" to be contained within the 
vertical dimension of the duct, the sound will propagate 
foorly. Frequencies corresponding to those longer wavel- 
engths will propagate poorly, if at all; higher frequencies 
With shorter wavelengths will fropagate much better. Ncrmal 
mode theory provides a tool for understanding these frofaga- 
tion effects [Ref. 3}. The transition from poor’ to good 
propagation is marked by the cutoff frequency. This 
frequency does not mark an abrupt end to propagation below 
which sound cannot be transmitted at all; rather it isa 
lower boundary below which propagation is less effective. 

Acoustic propagation in a sound channel is exactly anal- 
ogous to the propagation of radio waves in a ground-based 
duct . Urick [Ref. 2} provides an eguation from the theory 
of radio propagation for the maximum wavelength which will 
"fit" in such a duct. Dosso and Chapman [Ref. 1] use this 
equation to calculate the cutoff frequency for their 
80-meter thick shallow sound channel as about 250 Hz. 

Urick jives another eguation derived from radio frofaga- 


tion theory to calculate the cutoff frequency: [Ref. 2] 


=a" 
er =4 °/x IO H i oe 


where ae is the Maximum wavelength of the trapped sound 
in feet, and H isthe thickness of the duct in feet. 
Converted into metric units and using a sound velocity of 
1483 mys, this formula yields a cutoff frequency of about 


243 Hz for Dosso and Chapman's 80-meter channel. A third 
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formula iS given by Kinsler, Frey, Coppens and Sanders 
Pet. on 


2 
F ~ 0 x IO 3 
CO | ) /9 (egn 2.2) 


fOr [ ) in meters. This yields a cutoff frequency of about 
280 Hz for this channel. All of these calculations agree 
within useful ranges, Since the cutoff frequency does not 


mark an abrupt end to propagation. 


C. TRANSMISSION LOSS 


Several mechanisms act to decrease the energy of an 
acoustic wave aS the wavefront moves farther away from its 
source. The most obvious loss mechanism is spreading. Once 
an outwardly propagating wavefront reaches a reasonable 
distance from a sufficiently small source, the wavefront can 
be considered for all practical purposes to be spherical. 
Thus, the energy in the wave is spread out over an area that 
increases as the square of the range from the source. In 
decibels, spherical spreading alone would produce a trans- 


mission loss (TL) of 


T Lespuee CO lOO), R (eqn 2.3) 


If the sound is scmehow contained so that spreading may 
only take place horizontally, as is the case when the sound 
source is inside a channel or duct, the spreading is cylin- 


drical rather than spherical and the loss is yiven by 


TLey = 10 log, R (ean 204 


Other mechanisms cause additional energy to be lost. 


These attenuation mechanisms include scattering and 
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absorption. The degree of attenuation is strongly dependent 
on the frequency of the’ sound, With higher frequencies 
experiencing greater attenuation. Thorp [Ref. 4] develored 
extensive tables to describe the attenuation of sound as a 
function of freguency. 

Scund can be scattered from rough boundaries at the sea 
surface or bottom, or from discontinuities at the edges of 
the layers. Bubbles, marine organisms and suspended partic- 
ulate matter will alsc scatter sound [Ref. 3]. Scattering 
is generally more important for higher frequencies than for 
lower ones. 

Sound energy can also be lost through absorption. 
Viscosity effects in the interaction of water molecules as 
the sound wave passes account for one absorption mechanisn. 
Chemical reactions are also important sources of absorptive 
loss. Of the many chemical constituents of seawater, 
certain compounds appear to undergo continuous dissociation 
into ions and reassociation. This relaxation mechanism 
takes upd acoustic energy. The most important absorption 
mechanism for frequencies near 1000 Hz appears to be the 
borate-boric acid relaxation mechanism [Ref. 5]. 

In general, attenuation increases logarithmically with 
frequency. The attenuation losses due to such mechanisms 
are so small at low frequencies that the measurements must 
be made over extremely long ranges to detect any losses. 
The precision with which the measurements must’ be made and 
the inhomogeneities in the ocean further complicate the 
task. Urick [Ref. 2] provides an excellent discussion of 
absorption and attenuation. 

Two factors are operating with opposing effects on sound 
being prepagated through sound channels. First, attenuation 
mechanisms tend to cause losses of acoustic energy. Second, 
the focusing effect of the channel tends to reduce the 


losses. In the ideal case, a sound channel bounded above 
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and below by smooth interfaces between layers with no 
absorption or leakace should produce losses at about the 
rate due to-—cylindrical spreading alone. The mininun 
frequency which will be trapped well in achannel is the 
cutoff frequency. At frequencies near 1000 Hz, absorption 
mechanisms become impfertant, and scattering losses also come 
atom pl ay. Therefore it should ke possible to find scme 
optimum frequency for transmission of sound in a channel of 
given thickness, lying somewhere above the cutoff frequency 
and below the freguency range where absorption and. scat- 


tering losses become large. 
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A. DATA COLLECTION 


Desso and Chapman [Ref. 1] conducted an experiment in 
sound propagation within and below the shallow sound 
channel. During summer conditions in the northeast Pacific 
Ccean off the west coast of Canada, CFAV ENDEAVOUR took up a 
station and monitored data collection equipment. A single 
calibrated hydrophone was suspended below the vessel ata 
depth of 101 meters, approximately the shallow sound channel 
axis. An array of hydrophones was suspended at 417 meters 
near the deep sound channel axis. A second vessel, CVAF 
PARIZEAU, opened range from CFAV ENDEAVOUR, dropping explo- 
sive charges every 1.8 km (approximately one nautical mile). 
The charges were set to detonate at 98 meters, near the 
Shallcw sound channel axis. 

CKFAV PARIZEAU also measured sound velocity profiles at 
16-kilometer intervals over the track (Figures 3.1 through 
ees) - Later analysis revealed that the sound channel deep- 
ened over the portion of the track between eight and twenty 
kilometers, so that the sources were not located within the 
SSC and sound trapping was greatly reduced for sources deto- 
hated between those ranges. Section D below discusses steps 
that Dosso and Chapman took to account for this deepening of 
tae oSC. 

The received acoustic energy was measured as 1/3 octave 
-and averages for each shot for both the shallow and deep 
hydrorhones. Dosso and Chapman were able to distinguish 
bottom—interacting propagation paths by the different 
arrival times for shcts in the first 35 kilometers of the 


track, and these arrivals were subtracted from the resulting 


ZS 


energy levels. Bottom interactions could not be isolated 
from more distant signals. Propagation loss was calculated 
aS knewn source level minus the received level measured at 


the hydrofhones. 


Be. RESULTS AND DISCUSSION 


For the shallow sound channel [Ref. 1], measured propa- 
gation loss was plotted as a function of range for the 1/3 
octave kLands. The propagation loss demonstrated strong 
dependence on range, as expected. An optimum freguency for 
propagation of about 800 Hz was observed. 

Figure 3.4 shows the propagation loss measured by Dosso 
and Chapman for both source and receiver located within the 
Shallecw sound channel. The solid curve represents geone- 


trical spreading loss, calculated fron 


H = |O 109, (R- R.) (eqn 3.1) 


where range ne and ocean depth R are measured in meters. 
Figure 3.4 (a) represents propagation loss below the optinun 
frequency cf 800 Hz; Figure 3.4 (b) represents propagation 
loss above 800 Hz. The close agreement between the geomne- 
trical loss curve and the 800 Hz data over much of the range 
is noteworthy. 

Figure 3.5 displays the propagation loss within the 
Shallow sound channel for the optimum propagation frequency 
of 800 Hz, with the effects of chemical absorption removed. 
Again the solid line represents jeometrical spreading for 
comparison. Note the abrupt increase in losseS over the 
8-20 km range. Over this range the sound channel axis deep- 
ened, so that the sources were no longer within the duct. 

Figure 3.6 presents the same data plotted as a function 
of frequency for three-point range averages centered at 24 
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and 33 km The solid circles represent actual measurements; 
the open circles represent values adjusted upward to remove 
the effects of chemical absorption. The dotted lines repre- 
sent the expected lcss from geometrical spreading alone at 
eug Hz. This display shows more clearly the optimun 
frequency of 800 Hz and how closely the losses for this 
frequency approximate the losses due to geometrical 
spreading alone. It was clear that frequencies above and 
below 800 Hz experienced greater losses. DELEractionsor 
"leakage" from the shallow sound channel accounted for the 
losses at lower frequencies. High freguencies exrerienced 


losses due to Chemical absorption mechanisms and scattering. 


C. CCMPARISON BETWEEN SSC AND DSC PROPAGATION 


The deep sound channel trapped lower frequencies (longer 
wavelengths) than the shallow sound channel because the DSC 
is much thicker. Figure 3.7 demonstrates the effects of 
receiver depth on fropagation loss. The sources were 
located at 98 meters, near the shallow sound channel axis. 
One receiver was also near the Shallow sound channel axis at 
101 meters. The other receiver was positioned at 417 
meters, ‘ear the deer sound channel axis. At 160 Hz sound 
leaked easily out of the shallow sound channel but was well 
trapped within the deep sound channel. At.315 0iz 9 | JUS 
above the cutoff frequency for the shallow sound channel, 
propagation in the channel improved. At 630 Hz sound was 
effectively trapped in the shallow sound channel, and propa- 
gation loss at the deep receiver was significantly greater 
than at the shallow receiver. 

Dosso and Chapman plotted relative gain between the two 


receivers according to the relationship 


Relative Gain <P st a (eqn 3.2) 


SHALLOW 
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Figure 3.8 displays the results for a four-point range 
average centered at 24 km. The cutoff frequency of 250 Hz 
is clearly demonstrated. Below 250 Hz the SSC fails to trap 
the sound effectively, giving larger transmission losses 
than the DSC and a négative relative gain. Above 250 Hz the 
SSC traps the sound tetter, producing smaller transmission 
losses than the DSC and a positive relative gain. Larger 
propagation losses inthe SSC below 25 Hz result from 


surface decoupling effects. 


D. MCDELING PROPAGATION LOSS IN THE SSC 


Gecnetrical spreading is a simple model, but it is not 
able to predict the frequency and range dependence that 
Dosso and Chapman observed. A split-step parabolic equation 
(PE) model was used to calculate propagation loss Lkaseéed on 
the sound velocity profiles measured at ranges of zero, 16 
and 33 km along the track. [ Ref. 1]. Attempts to model 
propagation using only these three profiles, however, were 
unsatisfactory. The 16 km profile indicated that the axis 
of the £SC had deepened, so that the sources were no lilcnger 
contained in the channel. This is reflected in the increased 
propagation loss over the 5-20 km range in Figure 3.5. The 
improved propagation at and beyond 20 km indicated to Dosso 


and Chapman that the axis of the SSC must have risen to 


about 100 meters again. For this reason they used _ the 
profile recorded at 233 km for the 20 km range also. Losso 
and Chapman also truncated their four sound velocity 


frofiles at a depth of 600 meters, introducing an acousti- 
cally transparent bottom at that depth {Ref. 7]. This step 
had the practical effects of limiting tne Fourier transform 
matrix to a manageable size and shortening computational 
time somewhat. It may also have nad some effect on the 


predicted losses. 


ie 


With these four sound velocity profiles Dosso and 
Chapman's model predicted losses with an acceptakle degree 
of accuracy at 25 Hz and 630 Hz (the practical lower and 
upper limits of their model). But near the cutoff frequency 
of 250 Hz, the predicted losses were not in such good agree- 
ment with observed losses. The model did, however, predict 
strong dependence on both frequency and range as Dcssc and 
Chapman observed. Figure 3.9 shows the model results 


compared to actual measurements. 


E. SUMMARY OF DOSSO AND CHAPMAN'S EXPERIMENT 


Dosso and Chapman point out that the shallow. sound 
channel behaved like an acoustic waveguide with optimum 
Fropagation at 800 Hz, with increased losses at lower 
frequencies due to diffraction and at higher frequencies due 
to absorption and scattering [ Ref. 1]. The propagation was 
found to be very sensitive to changes in the environment. 
The PE model that Dosso and Chapman used was able to account 
for these range-dependent effects with good agreement with 
experimental results at low and near-optimum frequencies. 
At intermediate frequencies the model showed range and 
frequency dependent effects, but the predicted losses were 
greater than experimental results indicated. Dosso and 
Chapman emphasize that more accurate modeling of prcpagation 
losses requires a much finer sampling of the environment; 
that is, sound velocity profiles at much closer intervals 


than 16 kilometers. 
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25 Hz (Dosso and Chapman). 
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Dosso and Chapman reported good results [Ref. 1] using 
their split-step Parabolic Eguation (PB) model to predict 
propagaticn loss in a range-dependent environment where both 
a shallow sound channel and a deep sound channel existed. 
The first step in this study was to attempt to reproduce 
those results using the PE model available at Naval 


Fostgraduate School Mcnterey. 


Aj. MODEL DESCRIPTION 


The split-step Parabolic Equation (PE) acoustic trans- 
mission loss model was chosen for study because it preduces 
a humerical solution for an approximation to the wave equa- 
tion and is generally considered to provide the most accu- 
rate tool for predicting transmission loss. It is often 
used as a benchmark against which to compare the performance 
of other transmission loss models. It can deal with range- 
dependent envircnments (multiple sound velocity profiles) 
and with some interaction of sound with the sea floor. Its 
drawbacks are its relatively low-frequency capability, 
inability to handle all bottom loss conditions, and lack of 
flexibility in dealing with boundary conditions at the sea 
surface (the model uses a flat pressure-release boundary 
which does not account for the effects of wind and waves). 
This PE model is implemented on the IBM 3033 computer at the 
Naval Postgraduate School (NPS) and is accessed through an 
interactive program developed and managed by the schcol's 
oceanography department. 

The FORTRAN source code for this PE model was written by 
the Acoustic Environmental Support Detachment (AESD) of the 


a, 


Naval Ocean Research and Development Activity (NORDA) during 
the 1970's. A complete description of the physics and algo- 
rithms for this model are given in [Ref.- 6]. The wave equa- 


tion 


Py * ‘KS ‘4 ae P(r 2) — () (eqn 4.1) 


is an elliptic eguation which describes acoustic pressure as 
a function of range r and depth z. It can be approximated 


by a parabolic partial differential equation of the forn 


Y= (AB)Y a 


where 
a 
As ae 
—P oS) ee (eqn 4.3) 
Ko {ez 
and 





B = > (n= - |) (eqn 4.4) 


This solution makes the assum “ions that the sound is 
propagating radially outward from a source in a cylindri- 
cally symmetric medium and that the range ris very large 
compared to the size of the source and the wavelength cf the 
sound: these are the cylindrical symmetry and the far-field 


assumptions. The source is further considered to be oscil- 
lating harmonically at a discrete freyjuency. The ccean 
Surface is considered to be a flat pressure-release 


boundary, andthe acoustic field vanishes in the bottom at 
the maximum depth of the Fourier transform [Ref. 6]. 

The parabolic equation accounts for diffractior and 
other full-wave effects; it also allows solution for range- 


dependent environments. The parabolic eguation can be 
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solved numerically using the Tappert-Hardin split-step 
Fourier algorithm tc yield a solution for the entire 
acoustic pressure field as a function of range and depth 
[Ref. 6}. 

The PE solution is best where the source and receiver 
are arranged so that the path between the two makes a small 
angle with the hcrizcnotal. At larger angles, the group and 
phase velocities of the computed solution become subject to 
errors and produce inaccurate representations in the inter- 
ference pattern. This affects the accuracy of the predicted 
transmission loss fcr a given range. Computing time 
increases quickly with frequency. The numerical soluticn is 
limited by the practical Fourier transform size, which means 
that the input envirenment governs the extent to which the 
model can handle higker frequencies. For practical purfoses 
the PE model should be used to predict propagation loss for 
frequencies below 200 Hz along waterborne or shallow botton- 
kounce paths [{Ref. 6}. 

The FE model beccmes inappropriate for situations where 
the bottcm slopes steeply and where large sound speed gradi- 
ents cccur [{Ref. 6]. 


Be. VALILATION AGAINST DOSSO AND CHAPMAN'S DATA 


This study is concerned with waterborne paths and 
Shallcew angles. We are not specifically interested in 
bottom interactions. Even though the frequencies in which 
we are interested may be as high as 800 Hz [Ref. 1], we are 
willing to use computational time to examine a full-wave 
solution. Especially important for comparison with Dosso 
and Charman's experimental data is the fact that the PE 
model can utilize a range-dependent environment. The input 
parameters described in Appendix A are derived from [Ref. 1] 
and [ Ref. 7}. The resulting information and displays show 


reasonakle agreement with Dosso and Chapman's work. 
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Ccmpare Figures 4.1 through 4.5 with Figure 3.4. The 
smooth, downward trending curve represents transmission loss 
due to geometric spreading alone, calculated from equation 
3.1. The measured 1/3 octave rand levels of Figure 3.4 show 
Signs of smoothing because many frequencies have been added 
together, lessening the multipath interference effects that 
kecome apparent for discrete frequencies [Ref. 1}. 

It is apparent, however, that frequencies below the 
cutoff frequency of 250 Hz propagate poorly within the 
Shallow sound channel and that higher frequencies propagate 
better. An interesting outcome here is that the 2000 Hz 
experiences less transmission loss than the experimentally 
determined optimum frequency of 800 Hz at ranges greater 
than 17.5 kn. Note also that the three highest frequencies 
all show sudden steer increases in transmission loss in the 
17-20 km range and that the levels never really recover. 
This is the range at which Dosso and Chapman noted a change 
in their sound speed rfrofiles so that their sources were no 
longer within the shallow sound channel. 

Dosso and Chapman observed an increase in transmission 
loss over roughly the same ranges (Figure 3.4), followed by 
ketter propagation at ranges beyond 20 kn. They were able 
to model these results successfully by adding another sound 
velocity profile at 20 km, using the same profile they had 
measured at 33 kp. With the same input parameters the NPS 
PE model yielded a much steeper drop and did not indicate a 
return to levels expected from the experimental data. 

Dosso and Chapman's Figure 3.5 compares to Figure 4.6. 
For Figure 3.5 Dosso and Chapman plotted measured data, with 
losses due to chemical absorption removed. This figure 
clearly shows the close correspondence between the 
measurements for the optimum frequency of 800 Hz and the 
losses predicted for geometrical spreading alone. The only 


anomalies arise over the 8-20 km range, as discussed above. 
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The NPS model allows the user to select whether attenuation 
will be considered. Figure 4.6 represents the NPS PE model 
prediction for 800 Hz without attenuation. Again the steep 
drop in transmission occurs between 15 and 20 km. A compar- 
ison with Figure 4.2 shows that the two curves are essen- 
tially similar, but attenuation accounts for a 1-2 dB 
greater loss in Figure 4.2. 

Dessc and Chapman's experimental data yielded good 
results when propagation loss was plotted as a £unction of 
freguency for 24 and 323 km (Figure 3.6). Recall that Dosso 
and Charman measured 1/3 octave band levels and then calcu- 
lated three-point range averages to obtain each curve. The 
NPS PE modé€éi produces transmission loss predicticns for 
discrete frequencies. To obtain transmission loss for the 
desired ranges, values were read from the computer-generated 
PE model graphic output at the appropriate ranges and at one 
kilometer on either side. The three values for each range 
were averaged to produce the values used in plotting Figure 
Ce ae Again multipath interference effects can account for 
the fluctuations with range and frequency. While an optinun 
frequency of 800 Hz is easy to read from Figure 3.6, it is 
more difficult to state an optimum freguency from Figure 
4o7. Frequencies in the 500-700 Hz range seem to have 
propagated best. 

Dosso and Chapman compared results from shallow sound 
Channel and deep sound channel receivers in Figure 3.7. 
Figures 4.8 through 4.10 present the corresponding results 
from the NPS PE model. As expected for a frequency well 
above cutoff, the 630 Hz curves (Figure 4.8) show better 
propagation where both source and receiver were within the 
Shallow sound channel, but much greater loss (up to 10-20 
dB) where the receiver was in the deep sound channel. The 
two curves are in closer agreement for 315 Hz (Figure 4.9) 


as would be expected for a freguency near cutoff, where 
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trapping and leaking are better balanced. There is stilla 
10 dB rEetter propagation for the 15-17.5 km range when both 
source and receiver are in the channel. Figure 4.10 demon- 
strates significant transmission loss for 160 HZ, which is 
well telow cutoff. In this case the deep receiver experi- 
ences Ltetter propagation by 5-10 dB over almost the entire 
range, indicating that the shallow sound channel failed to 
trap this frequency effectively. 

It should especially be noted in Figures 4.8 and 4.9 
that ancmalous propagation is occurring in the 17-25 kn 
range. This may be due to the change in depth of the 
Shalicw sound channel that Dosso and Chapman noted and tried 
to correct by adding the sound velocity profile at 20 km. 
This ccmpensation was apparently much less successful with 
the NPS FE model than with Dosso and Chapman's model. 

Figure 4.11 corresponds to Figure 3.8. Relative gain is 
calculated according to equation 3.2, following Dosso and 
Chapman's procedure. Their results were good for 24 kn, 
Showing a sudden improvement above the 250 Az cutoff 
frequency when both source and receiver were in the SSC. 
The curve for the NPS PE model at the 24 km range does not 
produce good results, probably due to the failure of the 
added profile at 20 km to compensate for the changing depth 


Of the SSG. The 15 km curve, however, is taken before the 
SSC deepened and indicates better results. The 250 %8z 
cutoff frequency is clearly indicated. The NPS results, 


however, show unexpected peaks at 40-50 Hz and at 100 Hz, 
where the shallow sound channel receiver should not ke doing 
as well as the deer receiver. Other frequencies’ lelow 
cutoff behave as expected. 

Together Figures 4.12 through 4.15 correspond to Figure 
S63) 5 where Dosso and Chapman plotted their PE _ sodel 
predicticns alony with experimentally measured data. Note 


that Dosso and Chapman's 630 Hz curve showS a minimum at 
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about 18 kag, followed by better propagation at longer 
ranges. Their predictions correspond rather well to the 
measurements.— On the other hand, Figure 4.12 predicts 
several dB better propagation for 630 Hz than was observed 
for ranges less than 15 km and several dB worse than 
observed at ranges beyond 20 ka. Again the compensation 
using the 20 km profile was less successful for the NPS PE 
model than for Dosso and Chapman's model. The steer drop 
and failure to recover in the transmission loss curves are 
even more apparent in Figure 4.13 for 400 Hz and Figure 4.14 
fore 250 Hz. In Figure 4.15 for 25 Hz the observed values 
match guite well with predicted values out to about 15 kn. 
At that foint the prediction goes wrong again and predicted 


values fall 10-15 dB Eelow measured values. 


C. DISCUSSION OF NPS PE RESOLTS 


Discounting the multipath effects on transmission loss 
which produce extreme fluctuations in transmission loss with 
both range and frequency, it appears valid to use the NPS 


Parabolic Equation mcdel to estimate optimum freguency for 


sound propagation in a shallow sound channel. The plot of 
transmission loss versus frequency (Figure 4.7) Shows an 
optimum frequency range of 500-700 Hz. This is acceptable 


in ccmparison with LCosso and Chapman's measured optinun 
frequency of 800 Hz for the same set of sound velocity 
profiles. This is also good enough for active and passive 
sonar operators to use in selecting the best frequency 
ranges for their equipment. 

The NPS PE model was not, however, a good eStimatcr for 
absolute levels of transmission loss. The NPS PE curves 
Showed large differences in transmission loss levels when 
compared to Dosso and Chapman's measurements. For relative 


gain between deep and shallow receivers the NPS PE model 
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demonstrated an appropriate cutoff frequency around 250 Hz, 
Eut with unexpected peaks at frequencies well below cutoff. 

Scme of these anomalous results from the NPS PE model 
May arise from using full-derth sound velocity profiles. 
Reasoning from ray~tracing theory, a fully absorbing botton 
placed at 2500-2900 meters may have allowed sufficient depth 
for some deep propagation paths to be refracted back toward 
the receiver. The fully absorbing bottom at 600 meters 
which Dosso and Chapman used in their modeling may have 
eliminated propagation by these paths. Interference effects 
between these paths and the ducted propagation may account 
for scme of the differences between Dosso and Chapman's 
results and the NPS FE results. 
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Figure 4.2 Transmission Loss ws Range, 800 Hz (NPS PE Model). 
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Figure 4.4 Transmission Loss vs Range, 100 Hz (NPS PE Model). 


50 


REPRODUCED AT GOVERNMENT EXPENS 


REISE 


=f) 





Dihes ry ae ae one ay ree 
le, ‘OPA v3 wd od | Sobek ORG! 
__ _— TF Koy 8 ma eee Ce ? 
} : — : 5 mo he | 
LoT) Sow | Bice a Co eeke ete dws sca. 


Figure 4.5 Transmission Loss vs Range, 25 Hz (NPS PE Model). 
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Figure 4.6 Transmission Loss for 800 Hz, Absorption Removed 
(NPS PE Model). 
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Figure 4.7 Transmission Loss ws Frequency (NPS PE Model). 
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Pique 35 Shallow vs Deep Receiver, 630 Hz (NPS PE Model). 
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Figure 4.9 Shallow vs Deep Receiver, 315 Hz (NPS PE Model). 
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Figure 4.10 Shallow vs Deep Receiver, 160 HZ (NPS PE Model). 
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Figure 4.11 Relative Gain between Deep and Shallow Receivers 
(NPS PE Model). 
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Figure 4.12 Transmission Loss wsS Range, 630 Hz 
(NPS PE Model). : 
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Figure 4.13 Transmission Loss vs Range, 400 Hz 
(NPS PE Model). 
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Figure 4.14 Transmission Loss vs Range, 250 Hz 
(NPS PE Model). 
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Figure 4.15 Transmission Loss ws Range, 25 Hz 
(NPS PE Model). 
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V. TWO HIGH-LATITUDE SHALLOW SOUND CHANNELS 


A. OPTIMUM FREQUENCIES IN THE SHALLOW SOUND CHANNELS 


Since the NES PE model yielded acceptable estimates of 
the optimum frequency for Dosso and Chapman's shallow sound 
channel, the next step in this study was to apply the model 
to scme other Shallow sound channels. Appendix A contains 
the NFS PE model input information for two sound velocity 
profiles taken ina deep, high-latitude ocean basin during 
summer conditions (10 and 28 June). Figures 5.1 and 5.2 
show the SVPs. Note the presence of a surface layer and the 
absence cf a deep sound channel. 

The SSC for 10 June is 110 meters’ thick. Using equa- 
tions 2.1 and 2.2 we would expect cutoff frequencies of 148 
and 173 HZ for this channel. The SSC for 28 June is 105 
meters thick, yieiding cutoff frequencies of 159 and 186 Hz. 

A plot of transmission loss versus freguency is the best 
tool for determining optimum frequency. For these two 
profiles the source was modeled at 98 meters and the 
receiver at 100 meters depth, providing for propagation 
within the shallow scund channel. Agaix the transmission 
loss levels at each range were calculated from the NES PE 
model output graphics, by reading levels at the center 
ranges and one kilometer on either Side, then computing a 
straight arithmetic average of the three readings. In 
Figures 5.3 and 5.4 it is much more difficult to identify an 
optimum frequency for the two shallow sound channels. 

For the 28 June case (Figure 5.3) the transmission 
losses improve gradually and show optimum frequencies in the 
600-900 Hz range, but no readily identifiable effect in the 
vicinity of the cutoff frequency. Propagation at 
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frequencies below cutcff is much better than would have been 
expected from Dosso and Chapman's’ results. Figures 5.5 
through 5.8 are the NPS PE model results for several 
frequencies in the optimum range. Note that there is little 
difference in the overall transmission loss levels. 

The 10 June case (Figure 5.4) is even more difficult to 
analyze. Again there is no clear indication of cutoff 
frequency in the grarh. Note also that the curves for 15 
and 24 km are virtually flat for frequencies below 200 Hz. 
The two curves show that transmission loss is highly depen- 
dent upon freguency, with improved transmission at 315-400 
Hz and at several freguencies in the 600-1100 Hz range. 
Surprisingly, the transmission loss for 800 Hz on the 24 kn 
curve is some 3 dB greater than for 315 Hz. Figures 5.9 and 
5.10 are the NPS PE model results for the frequencies near 
optimum, and again there is not much apparent difference in 
the overall transmission loss levels. 

In an attempt to identify the optimum frequency for the 
10 June case more clearly, the investigator then introduced 
an artificially shallow bottom in hopes that this would 
Simplify the propagation paths. The 10 June sound velocity 
profile was truncated at the lower boundary of the shallow 
sound channel, with a fully absorbing bottom placed at 150 
meters. Source and receiver remained at 98 meters and 100 
meters respectively; no other input parameters were changed. 


Figure 5.11 displays the results. There is still no indica- 


tion of cutoff frequency, but two of the curves.’ dc show 
gradual igprovement in transmission with increasing 
frequency. The 15 and 24 km curves show smallest transmis- 
Sion losses in the 700-900 Hz range. Unfortunately the 33 


km curve shows smallest transmission loss at 600 Hz and 
large loss at 700 Hz. Figures 5.12 through 5.14 are NPS PE 
model results for frequencies near optimun. USI gmaierdr G1 = 


ficially introduced shallow bottom appears to have improved 
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the optimum frequency prediction somewhat for the 10 June 
case. 

Putting all three cases together and looking for a 
trend, we can tentatively identify a range of optimum 
frequencies around 600-900 Hz. The results for these 
Shallow sound channels are not as clear or conclusive as 
Dosso and Chapman's results, but they should be adequate for 


the fleet scnar operator to use in Setting up search plans. 


Be. OPTINUA FREQUENCY FOR PROPAGATION ACROSS THE LAYER 


The second part of this investigation involved the same 
sound velocity profiles witao the source still located at 98 
meters near the SSC axis, but with the receiver modeled at 
zO meters, well within the surface layer. Figures Same 
through 5.17 show the results. 

Figure 5.15 for 28 June shows peaks at 125 and 600 Hz in 
the transmission loss versus frequency curves for 15, 24 and 
33.” keg: Again the concept of cutoff freguency dces not 
appear to have much meaning in these cases. Figure 5.16 
Shows peaks at 160-200 Hz and at 700 Hz. For the truncated 
10 June profile, Figure 5.17 shows best propagation at 
75> V00-Eze There is a peak at 800 Hz, but the transmission 
loss is greater than at 75-100 Hz. 

Figures 5.18 thrceugh 5.20 show relative gain for our 


three cases, calculated according to the eguation 


Relative Gain =i a ia 


Saal 
LAYER ssc (“ae 


These high-latitude summer profiles do not show a cutoff 
frequency between the SC, and the surface layer. 
Propagation is always better when both source and receiver 
are contained within the shallow sound channel than when the 
source is in the SSC and the receiver isin the surface 


layer. 
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High-latitude Sound Velocity Profile, 10 June. 


65 


DEPTH (M) 
ZOO 2000 1500 1000 500 


SOOO 


SislOle, 


AOOO 


Figimesss2 





REPRODUCED AT GOVERNMENT EXPE 


SOURCE a» Se RECEIVERS 


1460 1470 41480 1490 #441800 += 1510~Ss«1520 
SOUND VELOCITY (M/S) 


High-latitude Sound Velocity Profile, 28 June. 
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Figure 5.3 Transmission Loss vs Freguency, 28 June 


(NPS PE Model). 
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Transmission Loss vs Frequency, 10 June 
(NPS PE Model). 
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nequre 5.5 Transsission Loss vs Range, 28 June, 600 Hz 
(NPS PE Model). 
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Figure 5.6 Transmission Loss vs Range, 28 June, 700 Hz 
(NPS PE Model). 
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Figure 5.7 Transpission Loss _ vs Range, 28 June, 800 Hz 
(NPS PE Model). 
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Figure 5.8 Transzission Loss _vs Range, 28 June, 900 Hz 
(NPS PE Model). 
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Figure 5.9 Transgkission Loss vs Range, 10 June, 700 Hz 
(NPS PE Mod el). 
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Figure 5.10 Transmission Loss_vs Range, 10 June, 800 Hz 
(NPS PE Model). 
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Figure 5.11 Iransmission Loss vs Frequency, 
Truncated Profile (NPS PE Model). 
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Figure 5.12 Transmission Loss vs Ran 


Truncated Profile, 


700 Hz 
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Figure 5.13 
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Figure 5.14 Transmission Loss vs Range, 


Truncated Profile, 


900 Hz (NPS PE Model). 
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Figure 5.15 Transmission Loss between Surface Layer and SSC, 
8 June (NPS PE Model). 
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Figure 5.16 Transmission Loss between Surface Layer and SSC, 
June (NPS PE Model). 
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Figure 5.17 Transmission Loss between Surface Layer and SSC, 
Truncated Profile (NPS PE Model). 
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Fagure 5.16 Relative Gain between Surface Layer and SSC, 
28 June (NPS PE Model). 
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Figure 5.19 Relative Gain between Surface Layer and SSC, 
0 June (NPS PE Model). 
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Figure 5.20 Relative Gain between Surface Layer and SSC, 
Truncated Profile (NPS PE Model). 
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A. MODEL DESCRIPTION 


The Fast Asymptotic Coherent Transmission loss nmecdel, 
version 9H (FACT 9H), like the split-step Parabolic Fquation 
model, is accessed through the Naval Postgraduate School's 
IBM 3033 computer. FACT currently provides the acoustic 
predictions for the Integrated Command ASW Predicticn Systen 
(ICAPS) on board U. S. Navy aircraft carriers and ASW 
Operations Centers (ASWOCS) worldwide. It is also the kasis 
for a number of acoustic products routinely provided by the 
Fleet Numerical Oceanographic Center, Monterey, to fleet 
users. Since the FACT model and its products are so widely 
distributed in the fleet, it is of interest to know how well 
FACT performs for the shallow sound channel case. 

The FACT 9H model as installed at the Naval Postgraduate 
School was developed for the Naval Ocean Research and 
Development Activity (NORDA) [Ref. 8]. It uses ray tracing 
theory with special corrections to improve the treatment of 
caustics, which would produce calculations indicating exces- 
Sive acoustic energy at certain ranges if not corrected. 
Unlike the PE model, FACT 9H accepts only one sound velocity 
profile, located at range zero. It makes no allowance for 
bottom tepography changing over the length of a track. The 
user may indicate the absorption characteristics of the 
bottom, or specify a table of values for bottom loss versus 
grazing angle. 

In addition, a semi-empirical expression is used to 
account for propagation ina surface duct. In this module 
the user may indicate the wave height, unlike the NFS PE 


model which makes no provisions for boundary conditions at 


oy) 


the surface. AS a fray-tracing model, FACT Q9Hlis “Hee 
expected to deal with waveguide phenomena such as diffrac- 
tion and leakage as successfully as the PE model does. The 
surface duct module is intended to remedy this problem sonme- 
what and to improve the acoustic prediction in the case of 
surface ducts. No correction has yet been applied for 
Shallcw sound cnannels. 

The investigation of the FACT 9H model proceeded along 
essentially the same lines as that of the PE model. The 
FACT 9H model was first tested against the experimental 
results cf Dosso and Chapman. Then the modei was applied to 
the 10 June and 28 June sound velocity profiles. Input 


parameters are descriked in Appendix A. 


Be. VALITATION AGAINST DOSSO AND CHAPMAN'S DATA 


Since FACT 9H does not allow for a range-dependent envi- 
ronment or for wave effects like diffraction and leakage, it 
was not expected to perform very well for the shallow sound 
Channel cbserved by Llosso and Chapman. Again the first test 
applied to the FACT 9H output was a plot of transmission 
loss versus frequency (Figure 6.1). As for the PE nmecdel, 
the levels were obtained by taking three-point range aver- 
ages centered at 24 kn. Since only the zero range sound 
velocity profile was used, there was no need to be ccncerned 
about tke effects of the shallow sound channel axis deep- 
ening over the eight to twenty kilometer range. 

Surprisingly, the transmission loss when both source and 
receiver were within the SSC showed a gradual decrease and 
optimum values at 800-900 Hz, consistent with Dosso and 
Chapman's actual measurements. Over the entire frequency 
band tested the FACT 9H model predicted losses between 70 
and 78 dR, while Dosso and Chapman observed losses between 
about 80 and 110 dB. The transmission loss in the optinun 
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frequency range is scme 10 dB less than Dosso and Chapman's 
measurements, but at least the optimum occurred at’ the 
proper frequencies. The transmission loss curve for the 
deep sound channel receiver showed only a fairly consistent 
increaSe with higher frequencies. 

Figures 6.2 and 6.3 are the FACT 9H transmission loss 
curves for the three frequencies near optimun. The shallow 
receiver curves (Figure 6.2) show more fluctuations, frob- 
ably due to multipath arrivals of sound reflecting from the 
flat ccean surface. The deep receiver curves (Figure 6.3) 
are smoothed by the partial absorption of sound in the ocean 
botton. 

The plot of relative gain (Figure 6.4) reiterates the 
information in Figure 6.1 ina way that should indicate the 
presence of a cutoff freyuency. We already know that the 
channel should exhibit a cutoff freguency at about 250 Hz, 
but there is no evidence for this in Figure 6.4. The 
Shallow receiver always shows better propagation than the 
deep receiver. This result is not surprising in that the 


FACT 9H model is not designed to handle such wave effects. 


C. THE HIGH-LATITUDE SHALLOW SOUND CHANNELS 


With good predictions for optimum frequency in the Dosso 
and Chapman case, the investigation proceeded to the 10 June 
and 28 June high-latitude, deep ocean sound velocity 
profiles. For the 28 June case (Figure 6.5) transmission 
loss in the shallow scund channel did indeed show a gradual 
rise and fall. The best frequency, however, was in the 
Vicinity of 400 Hz, much lower than the 600-900 Hz predicted 
by the PE model (Chapter 5, Section A). As one might 
expect, transmission losses between the SSC and the surface 
layer are always significantly greater than when the source 


and receiver are both contained in the SSC. 
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As for the 10 June case, FACT 9H results showed a slight 
but steady increase in transmission loss over the entire 
frequency range (Figure 6.6). The best freguency was, in 
fact, the lowest tested: TSW Le The transmission loss 
curve for the receiver located in the surface layer is only 
a few dB greater than the SSC curve. Recall that the NPS PE 
model also failed to froduce satisfactory results for the 10 
June case. Truncating the sound velocity profile at 150 
meters with a fully absorbing bottom seemed to improve the 
PE model prediction somewhat, but Figure 6.7 does not 
present much of a difference for the FACT 9H model. 

Perhaps the 10 June case iS a "weaker" shallow sound 
channel for modeling purposes. Dosso and Chapman's sound 
velocity profiles shcw a velocity gradient between 0.1436/s 
and 0.2071/s between the bottom of the surface layer and the 
axis of the shallow sound channel. The 28 June sound 
velocity profile has a velocity gradient of 0.1871/s. The 
10 June sound velocity profile, however, has a velocity 
gradient of only 70. 07377s. There may be a critical value 
for the velocity gradient between the surface layer and the 
SSC axis that will prcduce good results in the models. 

For now, however, it does not appear that the FACT 9H 
model is successful in predicting optimum frequencies for 
propagation of sound in shallow sound channels. In view of 
the extensive usage cf the FACT model to provide acoustic 
transmission loss products to fleet operations, this problem 


warrants further attention. 
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Figure 6.1 Iransmission Loss vs Frequency, 
Dosso and Chapman (NPS FACT 9H Model). 
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Figure 6.2 Transmission Loss in the Shallow Sound Channel, 
Dosso and Chapman (NPS FACT 9H Model). 
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Figure 6.3 Transmission Loss between SSC and DSC, 
Dosso and Chapman (NPS FACT 9H Model). 
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Figure 6.4 
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Relative Gain, Dosso and Chapman 
(XPS FACT 9H Model). 
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Figure 6.5 Transmission Loss vs Prequency, 28 June 
(NPS FACT 9H Model). 
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Figure 6.6 Transmission Loss vs Freguency, 10 June 


(NPS FACT 9H Model). 
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Figure 6.7 Transmission Loss vs Frequency, 
Truncated Profile (NPS FACT 9H Model). 
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_ VII. CONCLUSIONS 


A. RESULTS “ABD  DESCUSSION 


Two acoustic transmission loss models were examined to 
determine how well they predicted the optimum freguency for 
propagation of sound in shallow sound channels. The split- 
step Farabolic Equaticn (PE) model performed satisfactorily, 
while the Fast Asymptotic Coherent Transmission (FACT 9H) 
model was less satisfactory. 

The FE model has two advantages: it allows for a range- 
dependent environment and it isa full-wave approximation 
for the wave equation. These two factors enable the PE 
model to deal with wave phenomena like diffractior and 
ducting with a degree of success. The PE model predicted 
optimum frequencies cf 600-900 Hz for the shallow sound 
channels under study, consistent with the 800 Hz actually 
measured by Dosso and Chapman. The model seemed to produce 
ketter results when the sound velocity gradient between the 
overlying surface layer and the shallow sound channel was 
pronounced. The 28 June case, where the gradient was 
0.13871/s, yielded 600-900 Hz for optimum frequency, but the 
10 June case, where the gradient was only 0.0737/s, did not 
produce an optimum frequency. The 10 June predicticn was 
improved somewhat by introducing a fully absorbing bcttcm at 
the lewer boundary of the shallow sound channel. 

The FACT 9H model has its advantages, too: it is effi 
cient in terms of ccmputational time anditis the Navy's 
standard acoustic propagation loss model at the time of this 
writing. As a fray-tracing model, it was not expected to 
handle situations involving surface ducts or shallow sound 


channels as successfully as the PE model, even with the 
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semi-empirical surface duct module built into this version 
of FACT. The Dosso and Chapman sound velocity frofiles 
resulted in a surprisingly good prediction of 800-900 Hz, 
which agreed well with the experimental results. The 28 
June case, however, predicted an optimum of about 400 Hz; 
and the 10 June case did not predict any reasonable optimun. 
Again the strength of the sound velocity gradient may 
explain these different results. 

Of our two choices, the PE model, despite its long 
computation time, appears to be the better candidate for 
predicting optimum frequency of propagation ina shallow 
sound channel. The FACT model, despite its widespread use in 
the fleet, is not a good predictor for optimum frequency. 
The predictions for brcth models seem to be better when there 
is a large sound velocity gradient between the surface layer 


and the axis of the shallow sound channel. 


Be. RECCMMENDATIONS AND AREAS FOR FURTHER STUDY 


The fleet user has no easy access to a Parakolic 
Equation model to assess optimum freguencies for real-time 
environmental conditions. While this study predicts an 
cptimum range of 600-900 Hz where both source and receiver 
are contained within 80 to 110 meter thick channels, it 
would be helpful to draw up a table for other source- 
receiver depth combinations and channels of varying 
thickness. 

Several areas merit further investigation. First, this 
study uses a fully aksorbing bottom throughout. This is 
unrealistic. It would be of interest to assess botton 
losses more accurately and determine what difference, if 
any, this makes in the prediction of optimum frequency. 
Another area for investigation is the relative strength of 


Shallow sound channels. A strong channel, where the sound 
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velocity gradient between the bottom of the surface layer 
and the axis of the shallow sound channel, appears to trap 
sound more strongly and to produce a clearer prediction of 
optimum frequency. It would be helpful to find out how much 
of a velocity gradient is required to produce an optinunm 
frequency from each of the models, and whether the predicted 
optimum varies with the strength of the channel. 

Finally, other acoustic transmission loss models are 
available. The implicit finite difference (IFD) PE model is 
advertised as being able to handle boundary conditions at 
the ocean bottom more successfully than the split-step PE 
model. The IFD model should be investigated. The most 
recent revision of the FACT model, version 10A, is nearly 
ready fcr implementation. Technical documents [Ref. 9] 
indicate that FACT 10A will, among other things, handle 
ducting situations and boundary conditions better than FACT 
9H--this may affect performance for shallow sound channels. 
A third type of acoustic model is also available--the 
RAYMCLTE model. This model uses a combination of normal mode 
and ray-tracing techniques and may produce better results 
than FACT when applied to surface ducts and sound channels. 
The RAYMODE model has recently been designated as the Navy's 
new Standard acoustic propagation loss model, and its 
performance should be investigated for as wide a variety of 


environmental conditicns as possible. 
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APPENDIX A 


INBUT PARAMETERS FOR THE NPS PE AND FACT 9H MODELS 


Table I describes the input options available to the 


user through the Naval Postgraduate 
Parabolic Equation (PE) and Fast 
Transmission loss (FACT 9H) models. 


are the sound velocity profiles used in 


oy 


School's split-step 
Asymptotic Coherent 
Tables II through VI 
this study. 


TABLE I 
Input Parameters for the NPS PE and FACT 9H Models 


PARAMETERS PE MODEL FACT 9H 
Input units English English 
or metric Or metric 
Source depths One One 
Receiver depths Maxinum 20 One 
Freyuencies One Maximum 6 
Maximum range Unlimited 999 kilometers 
Range step Minimum .07 no Selected by 
or selected by prograr 
progran 
Number of SVPs Maximum 50 One 
Depth levels Maximum 50 Maximun 50 
in SVEFs 
Bottom lcss 1) Fully absorbing 1) Bottom type 
0 (reflecting) 
Eo (absorking) 
2) User's loss vs 2) User's loss 
grazing angle table VS grazing angle 
table 
=) User's depth _vs 3) Internal loss 
attenuation table vS grazing angle 
table 
Variable Maxinum 100 ranges/ No 
bottom depth depths 
Conrecticn for Available Available 
Spherical earth 
COrrecticn yror Available Available 
for vclumne 
attenuation 
Vertical bean 1-33 degrees No 
size of source 
Special ee Not reguired Yes 
for surface duc 
Wave height No Yes 
Critical angle Not required 1) Preset value 


8.65 degrees 


2) User input 
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TABLE 
Dosso and Chapman's Sound VWelocity Profile at 0 Kilometers 
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See Figure 3.1. 
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TABLE III 
Dosso and Chapman's Sound Velocity Profile at 16 Kilometers 
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DEPTH (M) 
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DEPTH (FT) 
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TABLE V 
High-latitude Sound Velocity Profile for 10 June 
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See Figure 5.1. 


TABLE VI 
High-latitude Sound Velocity Profile for 28 June 
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